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Abstract 

We review recent findings that the universe on its largest scales shows hints of 
the violation of statistical isotropy, in particular alignment with the geometry and 
direction of motion of the solar system, and missing power at scales greater than 
60 degrees. We present the evidence, attempts to explain it using astrophysical, 
cosmological or instrumental mechanisms, and prospects for future understanding. 
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1 Introduction 

The cosmological principle states that the universe is homogeneous and isotropic 
on its largest scales. The principle, introduced at the beginning of any cosmol- 
ogy course, is a crucial ingredient in obtaining most important results in quan- 
titative cosmology. For example, assuming the cosmological principle, cosmic 
microwave background (CMB) temperature fluctuations in different directions 
on the sky can be averaged out, leading to accurate constraints on cosmolog- 
ical parameters that we have today. However, there is no fundamental reason 
why statistical isotropy must be obeyed by our universe. Therefore, testing 
the cosmological principle is one of the crucial goals of modern cosmology. 

Statistical isotropy has only begun to be precision tested recently, with the ad- 
vent of first large-scale maps of the cosmic microwave background anisotropy 
and galaxy surveys. Extraordinary full-sky maps produced by the Wilkinson 
Microwave Anisotropy Probe (WMAP) experiment, in particular, are revolu- 
tionizingour ability to test the isotropy of the universe on its largest scales 
JH : S EI) Stakes are set even higher with the recent discovery of dark en- 
ergy that makes the universe undergo accelerated expansion. It is known that 
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dark energy can affect the largest scales of the universe — for example, the 
clustering scale of dark energy may be about the horizon size today. Simi- 
larly, inflationary models can induce observable effects on the largest scales 
via either explicit or spontaneous violations of statistical isotropy. 



2 Multipole vectors 



Multipole vectors are a new basis that describes the CMB anisotropy (or 
more generally, any scalar function on the sky) and are particularly useful in 
performing tests of isotropy and alignments. CMB temperature is traditionally 
expressed in harmonic basis, using the spherical harmonics Ye m . Copi, Huterer 
& Starkman (2003; (jjj) have introduced an alternative representation in terms 
of unit vectors 

2} = 5> m F £m (M) ~ f[(v^ • e) , (1) 

m i=i 



where v^' 1 ' is the i th multipole vector of the £ th multipole. (In fact the right 
hand side contains terms with "angular momentum" £ — 2, £ — 4 etc.; these 
are subtracted by taking, the appropriate traceless symmetric combination 
as described in (Ja) and (jfj).) In more technical language, Eq. (1) states the 
equivalence between a symmetric, traceless tensor of rank £ (middle term) 
and the outer product of £ unit vectors (last term). Note that the signs of all 
vectors can be absorbed into the sign of A^>. This representation is unique, 
and the right-hand side contains the familiar 2£ + 1 degrees of freedom - two 
dof for each vector, plus one for 

An efficient algorithm to compute the multipole vectors has been presented in 
(5) and is publicly available (7); other algorithms have been proposed as well 
( 8 ; 0; E3) • Interestingly, after the publication of the CHS paper (J^) , Weeks 
(|9|) pointed out that multipole vectors have actually first been invented by 



Maxwell ([lit ) more than 100 years ago! 



The relation between multipole vectors and the usual harmonic basis is very 
much the same as that between cartesian and spherical coordinates of stan- 
dard geometry: both are complete bases, but specific problems are much more 
easily solved in one basis than the other. In particular, we and others have 
found that multipole vectors are particularly well suited for tests of planarity 
of the CMB anisotropy pattern. Moreover, a number of interesting theoreti- 
cal results have been found; for example, Dennis (0) analytically computed 
the two-point correlation function of multipole vectors for a gaussian random, 
isotropic underlying field, while in Copi et al. (0) we have studied the relation 
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of multipole vectors to maximum angular momentum dispersion axes, max- 
ima/minima directions, and other related quantities that have been proposed 
to study the CMB. 



3 Alignments with the Solar System 

Armed with this new representation of the CMB anisotropy, we have set out 
to study the morphology of CMB anisotropies on large angular scales. Prior 
to our work, Tegmark et al. (flif ) found that the octopole is planar and that 
the quadrupole and octopole planes are aligned. In Schwarz et al. (jlih. we 
have investigated the quadrupole-octopole shape and orientation using the 
multipole vectors. Quadrupole defines two vectors and therefore one "oriented 
area" vector w^' 1 '" = x v^'. Octopole defines three multipole vectors 
and therefore three normals. Hence there are a total of four planes determined 
by the quadrupole and octopole. 

In Schwarz et al. we found that 

• the normals to these four planes are aligned with the direction of the cosmo- 
logical dipole and with the equinoxes at a level inconsistent with Gaussian 
random, statistically isotropic skies at 99.95% C.L.; 

• the quadrupole and octopole planes are orthogonal to the ecliptic at the 
98.5% C.L.; 

• the ecliptic threads between a hot and a cold spot of the combined quadrupole 
and octopole map, following a node line across about 1/3 of the sky and 
separating the three strong extrema from the three weak extrema of the 
map; this is unlikely at about the 95% C.L.; 

• the four area vectors of the quadrupole and octopole are mutually close (i.e. 
the quadrupole and octopole planes are aligned) at the 99.9% C.L. 

(These numbers refer to the TOH map; other maps give similar results as 
Table 3 of Ref. (0) shows.) While not all of these alignments are statistically 
independent, they are clearly surprising, highly statistically significant (at 
> 99.9% C.L.), and unexpected in the standard inflationary theory and the 
accepted cosmological model. 

Particularly puzzling are the alignments with the solar system features. CMB 
anisotropy should clearly not be correlated with our local habitat. While the 
observed correlations seem to hint that there is contamination by a foreground 
or perhaps scanning strategy of the telescope, closer inspection (see the next 
two Sections) reveals that there is no one obvious way to explain the observed 
correlations. 
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L=2+3 alignments 




Fig. 1. Quadrupole and octopole [l = 2 and 3) of the WMAP sky map in galac- 
tic coordinates, shown with the ecliptic plane, the supergalactic plane (SGP), the 
equinoxes and the cosmological dipole. We also show the four normals to the 
planes defined by vectors that describe the quadrupole and octopole temperature 
anisotropy; one normal is defined by the quadrupole and three by the octopole. Note 
that three out of four normals lie very close to both the equinoxes and the dipole 
direction. The probability of these alignments being accidental is about one part in 
a thousand. Moreover, the ecliptic plane traces out a locus of zero of the combined 
quadrupole and octopole over a broad swath of the sky — neatly separating a hot 
spot in the northern sky from a cold spot in the south. These apparent correlations 
with the solar system geometry are puzzling and currently unexplained. 



Our studies (see (p)) indicate that the observed alignments are equinoctic/ecliptic 
ones (and/or correlated with the dipole direction), and not alignments with 
the Galactic plane: the alignments of the quadrupole and octopole planes with 
the equinox/ecliptic/dipole directions are more significant than those for the 
Galactic plane. This conclusion is supported by the foreground analysis (see 
the next Section). Moreover, it is remarkably curious that it is precisely the 
ecliptic alignment that has been found on somewhat smaller scales using the 
power spectrum analyses of statistical isotropy (jloT : llfih. 

Finally, it is important to make sure that the results are unbiased due to 
unfairly chosen statistics. We have studied this issue extensively in (0). Two 
natural choices of statistics which define ordering relations on the three dot- 
products between the quadrupole and octopole area vectors A iy each lying in 
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the interval [0, 1], are: 



S=^(A 1 + A 2 + A 3 ), and 



T=l 



1 r 



[1-A 1 ) 2 + (1-A 2 ) 2 + (1-A ," 2 



Both S and T can be viewed as the suitably defined "distance" to the vertex 
(A 1 ,A 2 ,A 3 ) = (0,0,0). A third obvious choice, (A?+A|+A§)/3, is just 2S-T. 
To test alignment of the quadrupole and octopole planes (or associated area 
vectors) we quoted the S statistic numbers; T gives similar results. 

To test alignments of multipole planes, we define the plane as the one whose 
normal, h, has the largest dot product with the sum of the area vectors (0). 
Since uTj-n is defined only up to a sign — Wi is headless — we take the absolute 
value of each dot product. Therefore, we find n that maximizes 

<S = — ^ \wi • n\ , (3) 



where Ni is the total number of area vectors considered. Alternatively, gen- 
eralizing the definition in ffl3h . one can find the direction that maximizes the 
angular momentum and compare the maximal angular momentumffor the 
quadrupole plus octopole) with that from simulated isotropic skies (J6|) 



J2 e m =-e 777,2 




2 


L^m=—e 


dim 


2 



(4) 



which gives similar results as the S statistic for the alignment of £ = 2 and 
£ = 3. 



4 Foregrounds 



While the statistical significance of the observed vs. the Galactic alignments 
of the quadrupole and octopole, by itself, may not not sufficient to rule out 
the Galactic contamination, we have explored several lines of reasoning which 
suggest that Galactic foregrounds are not the cause of the alignments (see also 
studies by (0; El 03)). 



First, we have tried adding (or subtracting) known, measured Galactic con- 
tamination to WMAP maps and observing how the multipole vectors move 
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Galactic cut (degrees) Ecliptic cut (degrees) 

Fig. 2. Quadrupole-octopole probabilities for the TOH map for an increasingly larger 
isolatitude cut of ± (degrees shown), performed symmetrically around the Galactic 
plane (left panel) or the ecliptic plane (right panel). We consider the S statistic 
probabilities applied to the ecliptic plane alignment — the solid line is the mean 
value, while the dark and light regions represent 68% C.L. and 95% C.L. regions, 
respectively, from 1000 realizations of reconstructed a£ m coefficients. The dashed 
line denotes the probability obtained from the full-sky map, corresponding to the 
case of zero cut. Adopted from (0). 

. In the large- foreground limit, the quadrupole vectors move near the z- 
axis and the normal into the Galactic plane, while for the octopole all three 
normals become close to the Galactic disk at 90° from the Galactic center. 
Therefore, as expected Galactic foregrounds lead to Galactic, and not ecliptic, 
correlations of the quadrupole and octopole. 

Second, in (0), we have shown that the known Galactic foregrounds possess a 
multipole vector structure very different from that of the observed quadrupole 
and octopole. The quadrupole is nearly pure I22 m the frame where the z-axis 
is parallel to the dipole (or w^ 2 ' 1 ' 2 ^ or any nearly equivalent direction), while 
the octopole is dominantly Y33 in the same frame. Mechanisms which produce 
an alteration of the microwave signal from a relatively small patch of sky - 
and all of the recent proposals fall into this class — are most likely to produce 
aligned Y20 and Y 30 (essentially because the multipole vectors of the affected 
multipoles will all be parallel to each other, leading to a Yi in this frame). 

Most of the results discussed so far have been obtained using reconstructed 
full-sky maps of the WMAP observations (|l|; [k| Ecft . Results with the re- 
constructed full-sky map in the presence of the sky cut is shown in Fig. (2): 
even with a cut of a few degrees (iso- latitude, for simplicity), the errors in the 
reconstructed anisotropy pattern, and the directions of multipole vectors, are 
too large to allow drawing quantitative conclusions about the observed align- 
ments. Figure 2 does show, however, that the cut-sky alignment probabilities, 
while very uncertain, are consistent with the full-sky values. Ultimately, one 
will want to check for the \ow-£ alignments on Markov chain Monte Carlo 
maps, where realizations of the reconstructed the anisotropy pattern over the 
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whole sky are based on the observations outside of the Galactic cut. While in 
principle straightforward (see e.g. (1211)). the key issue in this approach that 
requires considerable care is modeling of the foregrounds. 



5 Quest for an explanation 



Understanding the origin of CMB anomalies is clearly important, as the ob- 
served alignments of power at large scales are inconsistent with predictions of 
standard cosmological theory. A number of authors have attempted to explain 
the observed quadrupole-octopole correlations in terms of a new foreground 
— for example the Rees-Sciama effect (Jiih. interstellar dust (24), local voids 
(jH),|26j). Most if not all of these proposals have a difficult time explaining 



the anomalies without severe fine tuning. For example, Vale (|27l ) cleverly sug- 
gested that the moving lens effect, with the Great Attractor as a source, might 
be responsible for the extra anisotropy; however Cooray & Seto (28) have ar- 
gued that the lensing effect is far too small and requires too large a mass of 
the Attractor. 



In Gordon et al. (|29l ) we have explored the alignment mechanisms in detail, 
and studied additive models where the temperature is added to the intrinsic 
temperature 

-^observed ( A) ^Intrinsic 

(n) + T add (n) (5) 



where T a dd(n) is the additive term - perhaps contamination by a foreground, 
perhaps an additive instrumental or cosmological effect. We have shown that 
additive modulations of the CMB sky that ameliorate the alignment problems 
tend to worsen the overall likelihood at large scales (they still may pick up 
positive likelihood contribution from higher multipoles). The intuitive reason 
for this is that there are two penalties incurred by the additive modulation. 
First, since the temperature at large scales is lower than expected, one typically 
needs to arrange for an accidental cancellation between T intrinsic and T add . 
Second, certain ag m in the dipole frame are observed to be suppressed relative 
to the expectation (see Table I in Ref. (|29h) - but zero is actually the most 
likely value of any given a,£ m , so likelihood with Ti ntr i ns i C is again penalized. 

Instead, the multiplicative mechanisms, where the intrinsic temperature is 
multiplied by a spatially varying modulation, are more promising. As a proof 
of principle, we suggested a toy-model modulation 

T b S erved(n) = / [1 + W 2 Y 20 (n)} 

-^Intrinsic ( n )> (6) 
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-18.2 ■ — 18.2(iK -33.7 — 33.7|xK 

quadrupole octopole 



Fig. 3. A realization of the multiplicative model where the quadrupole (left column) 
and octopole (right column) exhibit an alignment similar to WMAP. First row: 
intrinsic (unmodulated) sky from a Gaussian random isotropic realization. Second 
row (single column): the quadrupolar modulation oc — [1 — 7l2o(n)] in the dipole 
direction. Third row: the modulated sky of the observed CMB. Fourth row: WMAP 
full-sky quadrupole and octopole. 

(where the modulation is a pure Y 2 q along the dipole axis), we have shown that 
the likelihood of the WMAP data can be increased by a factor of exp (16/2) 
and, at the same time, the probability of obtaining a sky with more align- 
ment (e.g. higher angular momentum statistic) is increased 200 times, to 45%. 
(Spergel et al. 2006 (J4j) thereafter did a similar study, generalizing the multi- 
plicative modulation to eight free parameters corresponding to all components 
of the dipole and quadrupole and finding the highest likelihood fit; see their 



Finally we have considered a possibility of an imperfect instrument, where the 
instrumental response to the signal T(n) is nonlinear 



Since the biggest signal on the sky is the dipole (of order mK), leakage of about 
1% (i.e. a>i ~ «2 ~ 0.01), if judiciously chosen, can produce the quadrupole 
and octopole that are as observed and are aligned with the dipole. Unfortu- 
nately (or fortunately!), WMAP detectors are known to be linear to much 



Fig. (26)). 




(7) 
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better than 1%, so this particular realization of the instrumental explanation 
does not work. As an aside, note that this type of explanation needs to assure 
that the higher multipoles are not aligned with the dipole/ecliptic, and more- 
over, requires essentially no intrinsic power at large scales (that is, even less 
than what is observed). 



6 Missing angular power at large scales 



Spergel et al. have found that the two point correlation function, C(n-n') = 
(T(n)T(n')), nearly vanishes on scales greater than about 60 degrees, contrary 
to what the standard ACDM theory predicts, and in agreement with the same 
finding obtained from COBE data about a decade earlier (fjof). Using the 
statistic 

1/2 

S 1/2 = I [C(B)) 2 d{co S 9). (8) 
-l 



Spergel et al. found that only 0.15% of the elements in their Markov chain of 
ACDM model CMB skies had lower values of S±/2 than the observed sky. 

We have revisited the angular two point function in the 3-yr WMAP data 



in Ref. (|3l|). We found that the two-point function computed from the var- 
ious cut-sky maps shows an even stronger lack of power, now significant at 
the 0.03%-0.15% level depending on the map used; see Fig. (4). However, we 
also found that, while C{6) computed in pixel space over the unmasked sky 
agrees with the harmonic space calculation that uses the pseudo-Q estima- 
tor, it disagrees with the Ce obtained using the maximum likelihood estimator 
(advocated in the 3rd year WMAP release (0)). The MLE-based Cg lead to 
C{6) that is low (according to the S\/2 statistic) only at the 8% level. This is 
illustrated in the right panel of Fig. (4). We are concerned that the full-sky 
maximum-likelihood map making algorithm is inserting significant extra large 
angle power into precisely those portions of the sky where we have the least 
reliable information. Clearly, the definitive judgment of the large-angle power 
has not yet been made. 

Finally, here we note that the vanishing of power is much more apparent in 
real space (as in C(8)) than in multipole space (as in Cg). The harmonic- 
space quadrupole and octopole are only moderately low (e.g. (|2lh). and it is 
really a range of low multipoles that conspire to make up the vanishing C(9). 
Therefore, theoretical efforts to explain "low power on large scales" should 
focus to explain the low C{6) at 9 > 60 deg. 
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9 (degrees) 9 (degrees) 



Fig. 4. Two point angular correlation function, C{9) = T(e\)T(e2) e computed in 
pixel space, for three different bands masked with the kpO mask. Also shown is the 
correlation function for the ILC map with and without the mask, and the value 
expected for a statistically isotropic sky with best-fit ACDM cosmology together 
with 68% error bars. Left panel: year 1 results. Right panel: year 123 results. Even 
by eye, it is apparent that masked year 123 maps have C{6) that is consistent with 
zero at 9 > 60 deg, even more so than in year 1 maps (at the 0.03%-0.15% level 
depending on the map used). We also show the C(9) computed from the "official" 
published Ci, which (at i < 10) are the pseudo-C^ in year 1, the and MLE Ci in 
year 123. Clearly, the maximum likelihood estimator (MLE)-based Cg, as well as 
C{9) computed from the full-sky ILC maps, are in significant disagreement with the 
angular correlation function computed from cut-sky maps. 

7 Discussion and Future Prospects 



If indeed the observed t = 2 and 3 CMB fluctuations are not cosmological, one 
must reconsider all CMB results that rely on low is, such as the normalization, 
A, of the primordial fluctuations and any constraint on the running dn s /d\og k 
of the spectral index of scalar perturbations. Moreover, the CMB-galaxy cross- 
correlation, which has been used to provide evidence for the Integrated Sachs- 
Wolfe effect and hence the existence of dark energy, also gets contributions 
from the lowest multipoles (though the main contribution comes from slightly 
smaller scales, i ~ 10). Finally, it is quite possible that the underlying physical 
mechanism does not cut off abruptly at the octopole, but rather affects the 
higher multipoles. Indeed, several pieces of evidence have been presented for 
anomalies at I > 3 (j22[ flST) : see also relevant work in Refs. (0; EE EE EE EE 

E3E3). 

So far no convincing explanation has been offered. In fact, a no-go argument 
has been given by Gordon et al. (|29l). reasoning that additive mechanisms for 
adjusting the intrinsic CMB anisotropy lead to a lower likelihood at low i than 
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the observed sky. Therefore, it appears that a multiplicative mechanism is at 
work, whether it is astrophysical, instrumental or cosmological. 

While the further WMAP data (4-year, 8-year etc) is not expected to change 
any of the observed results, our understanding and analysis techniques are 
likely to improve. Much work remains to study the large-scale correlations 
using improved foreground treatment, accounting even for the subtle system- 
atics, and in particular studying the time-ordered data from the spacecraft. 
The Planck experiment will be of great importance, as it will provide maps of 
the largest scales obtained using a very different experimental approach than 
WMAP — measuring the absolute temperature rather than temperature dif- 
ferences. Polarization maps, when and if available at high enough signal-to- 
noise at large scales (which may not be soon), will be a fantastic independent 
test of the alignments, and in particular each explanation for the alignments, 
in principle, also predicts the statistics of the polarization pattern on the sky. 

The quest for an answer has whetted the appetite of cosmologists to under- 
stand the structure of the universe on its largest scales. 
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